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ABSTRACT 

In order to extract information about behavior of polarized gluons in the nu- 
cleon, charmed hadron productions, i.e. D* + meson and A+ baryon productions, 
are studied in polarized pp reactions at BNL-RHIC energy. For these processes, 
the spin correlation asymmetry Dll between the target proton and the produced 
charmed hadron, and its statistical sensitivity 5 Dll are calculated. From analyses 
on these processes, we found that the pseudo-rapidity distribution of D LL in the 
limited transverse momentum region is quite effective for distinguishing the model 
of polarized gluons as well as the model of spin-dependent fragmentation functions. 
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1. Introduction 

Though the proton is one of the most familiar nucleons, the origin of its spin is still 
an open question. About 15 years ago, the European Muon Collaboration (EMC) 
reported surprising data on the polarized structure function of proton gi(x, Q 2 ) [1], 
indicating that the spin of the proton cannot be described by the naive quark model. 
Disagreement of the result extracted from the EMC data with theoretical predic- 
tions has been called the "Proton Spin Puzzle" , which is not yet solved well in spite 
of many follow-up experiments and a great development of theoretical analyses. It 
is still one of the challenging themes in nuclear and particle physics [2]. Accord- 
ing to the parton-model concept inspired by quantum chromodynamics(QCD), the 
proton spin satisfies the following sum rule; 



are the spin carried by quarks, anti-quarks and gluons, respectively. In the right- 
hand side of Eq.(l), integration is performed over all Bjorken-x region. < L > q + g 
represents the orbital angular momentum of quarks and gluons in the proton. An 
important thing is to know how large each of these component is and then to 
understand, based on QCD, what the underlying dynamics of this sum rule is. 
These days the behavior of valence quarks in the proton has been considerably well 
understood with great efforts. However, the knowledge of gluons in the proton 
is still poor. Therefore, to solve the proton spin puzzle, it is very important to 
obtain a good information about gluon polarization in the proton. So far, the 
gluon distribution in the proton has been studied mainly through deep inelastic 
scattering (DIS) of polarized leptons (electrons and muons) off polarized nucleons. 
However, now we are in a new stage; the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory (BNL) which gives us another unique place 
for probing the proton structure via polarized proton-polarized proton scattering 
has started to work. For the RHIC experiment, the following processes are already 
proposed to study the behavior of polarized gluons in the proton [3]; 

• High-p T ("prompt") photon production, pp — > jX, 

• Jet production, pp — > jet(s)X, 




(1) 



where \ on the left-hand side is a spin of proton, while Aq(x), Aq(x) and Ag(x) 



• Heavy-flavor production, pp — ■> ccX, bbX, 

where the particle with an arrow means that it is polarized. 

In addition to these processes, we propose other interesting processes, i.e. the 
polarized charmed hadron productions in the polarized proton-unpolarized proton 
collision: 



by expecting that they also could be observed in the forthcoming RHIC experi- 
ment. In the lowest order of QCD, a charm quark being one of the constituents 
of charmed hadrons is produced via gluon-gluon fusion (gg — > cc) and quark-anti- 
quark annihilation (qq — > cc) processes in proton-proton collisions (Fig.l). Then, 
if we could separate kinematically a contribution of gg — > cc from qq — > cc and 
pick up gg — > cc clearly, we could extract a good information about gluon polariza- 
tion by observing the spin correlation asymmetries between the initial proton and 
these produced charmed hadrons, since the cross section of gluon fusion is directly 
proportional to the gluon distribution. 

Concerning the A+ production, in a couple of years ago we analyzed the spin 
correlation between the initial proton and produced A+ in the diffractive regions of 
the process pp — > pA+X, where Pomeron interactions play an important role [4] and 
pointed out that the process is promising for extracting the polarized gluon distri- 
bution. In this paper, we develop those analyses to another important kinematical 
region, i.e. the central region of collision. The present analysis is complementary to 
those previous analyses and should give us additional information on the polarized 
gluon distribution in the different kinematical region " 1 . 

This paper is organized as follows. In the next section, we explain why we focus 
on charmed hadron production in this work. Then, we introduce spin correlation 
observable and present a theoretical formulation of our calculation in Sec. 3. In 
section 4, we show the result of our numerical calculation. Finally, Sec. 5 is devoted 
to conclusion, including discussion and summary. 

111 The preliminary analysis of the A+ production was done in rcf [5]. 
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Figure 1: Diagrams for pp — > htX at the lowest order (left) and its subprocesses 
(right), where h+ means a polarized charmed hadron with positive charge. Dashed- 
lines in the left figure stand for partons which are emitted from initial protons. 

2. Polarized charmed hadron production in pp 
collisions 

As described above, we will focus on two processes of charmed hadron production; 
one is polarized D* + meson production and another is polarized A+ production. 
Then, what is so special about these processes? and why are they so interesting? 
This is because in these processes, we can regard the spin of the produced charmed 
hadron to be same as the one of the charm quark as described below; 
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Since the D* + meson is in the 3 Si state in the non-relativistic quark model, 
the spin of D* + meson is carried by a charm quark and an anti-down quark 
whose spins are combined in parallel. Therefore, the spin direction of the D* + 
meson and the spin direction of the constituent charm quark are expected to 
be same. 

• baryon 

The A+ is composed of a heavy charm quark and antisymmetrically combined 
light up and down quarks. Hence, the A+ spin is basically carried by a charm 
quark. 

Furthermore, it is expected to be very rare for a produced charm quark to change 
its spin direction during its decay into a charmed hadron, since the charm quark 
is significantly heavy and the spin flip interaction being proportional to l/m c (m c 
is charm quark mass) is small. Therefore, if these processes are originated from 



gluon fusion, spin-dependent observables for produced charmed hadrons directly 
depend on the polarization of gluons in the initial proton and thus, observation of 
the polarization of the produced charmed hadron could give us a good information 
about polarized gluons in the proton. 

Concerning the above discussion, it is worth while to comment on contribu- 
tion of light-quark fragmentation to charmed hadron production. In the case of 
charmed hadron production by light-quark fragmentation, a light-quark produced 
in a hard parton subprocess must pick up a charm quark created from vacuum to 
make a charmed hadron. However, the probability of charm-quark pair creation 
is extremely small, estimated to be 1CT 11 compared to the one for uu and dd pair 
creation [6]. Therefore, in charmed hadron production we can safely neglect the 
contribution of light-quark fragmentation. This is rather different from the case of 
A production [7] , where the probability of ss pair creation from vacuum cannot be 
neglected. 



3. Spin correlation observables and cross sections 

To study the polarized gluon distribution in the proton, we introduce the spin cor- 
relation asymmetry between the polarized target-proton and the produced charmed 
hadron defined by; 

da ++ — da + _ + da — da_ + 
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da ++ + d<7 + _ + da + da_ + 

dAa/dX 



da/dX ' 



(X = p T or 77), (4) 



where da + -, for example, denotes the spin-dependent differential cross section with 
the positive helicity of the target proton and the negative helicity of the produced 
charmed hadron. p T and 77, which are represented by X in Eq.(4), are transverse 
momentum and pseudo-rapidity of the produced charmed hadron, respectively. 

In order to perform our analyses, we take the proton-proton center of mass 
system(CMS). In this system, four momentum p x of each particle % is defined as 



follows; 



Js / 4m 2 

P A/B = ^(^ TP: 0) , 0=^1-—* 

P hc = (E P L , P T ) ( 5 ) 



{^m 2 hc +p T 2 cosh 2 r], p T sinhri, p T ), 



Pa/b X a/bPA/B' Pc 5 

where the parameters in parentheses denote the energy, the longitudinal momen- 
tum and the transverse momentum, in this order, stands for the mass of the 
particle i (i = p for proton and % = h c for A+ or D* + ). x a / b and z are momentum 
fractions of the parton to the proton and of the charmed hadron to the charm 
quark, respectively. Here, notice that we define the momentum of unpolarized and 
polarized proton as p A and p B , respectively. In addition, we regard the direction 
of p B as the positive z-axis. 

According to the quark-parton model, p T - or ^-dependent differential cross sec- 
tion, dAa/dX, is given by, ^ 2 

dAa /" ymax f 1 f 1 T n2 , A , ^dAa(gg^cc) 
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(6) 



dt 

x AD nt/ .(z)Jdx a dx b dY, (X, Y = p T or r\ [X ^ Y)) 

where g ga /p A (x a ,Q 2 ) and Ag$ (x b , Q 2 ) are the unpolarized and polarized gluon 
distribution functions, respectively, and q qa (q b )/ PA (p B )(x a{b) , Q 2 ) and Aq qa(<ib)/pA{l}B) (x a{b) , Q 2 ) 
denote the unpolarized and polarized distributions, respectively, of the quark and 
anti-quark. AD^+ ^(z) represents the spin-dependent fragmentation function of 
the outgoing charm quark decaying into a polarized charmed hadron. The spin- 
dependent differential cross sections of the subprocess is given by 113 
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* 2 rn Eq.(6), X or Y mean p T or r\. Thus, if X is p T , then Y is 77, and vice versa. 
" 3 These algebraic calculations are carried out using FORM [8]. 
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and the Jacobian, J, which transforms the variables z and t into p T and 77, is given 
by 
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zs^Jm 2 + p^cosh 3 ?7' 
where we defined the following variables; 

S = X a XbS, t\ = — t, U\ = U, 

z z 

with 
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In Eq.(6), the minima of a; a and Xb are given by 
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The unpolarized differential cross section for this process is calculated by replac- 
ing spin-dependent functions Ag(x), Aq(x), AD(z) and by spin-independent 
ones g(x), q(x), D(z) and ^, respectively " 4 and we use it to estimate the denom- 
inator of Dll. 



4. Numerical calculation 

To carry out a numerical calculation of D LL , we used, as input parameters, 
to c = 1.20 GeV, Trip = 0.938 GeV, m D , + =2.01GeV and m A + = 2.28 GeV [10]. 
In addition, we used the AAC [11] and GRSV01 [12] parameterization models for 
the polarized parton distribution function and the GRV98 [13] model for the un- 
polarized one, and set the scaling variable Q 2 as Q 2 = p 2 T - Though both AAC 

" 4 The spin-independent subprocess differential cross section ^ was given in Ref. [9] 



and GRSV01 models excellently reproduce the experimental data on the polarized 
structure function of nucleons gi(x), the polarized gluon distributions for those 
models are quite different. In other words, the data on polarized structure function 
of nucleons g±(x) alone are not enough to distinguish the model of gluon distribu- 
tions. Since the process is semi-inclusive, it is necessary to know the fragmentation 
function of a charm quark to D* + or A+ to carry out numerical calculations. For 
the spin-independent fragmentation function Dh c / C (z), we use the Peterson frag- 
mentation function which was proposed by Peterson et al. [14] a long time ago and 
has been widely used for phenomenological analyses. This fragmentation function 
include one free parameter which is determined by experimental analysis. In this 
work, we take e£)»=0.078 [15] and eA c =0.25 [16] for D* + production and A+ pro- 
duction, respectively. Furthermore, for the spin-dependent fragmentation function 
of a charm quark decaying into polarized charmed hadrons, AD^^(z) , we have no 
information about it at present, since we have no data on polarized charmed hadron 
productions. Therefore, here we simply assume the spin-dependent fragmentation 
function to be parametrized [17] as 



where Ch c / C is a scale-independent spin-transfer coefficient. We consider the fol- 
lowing two models: 

(A) Ch c / C = 1 (non-relativistic quark model) 

(B) Ch c / C = z (Jet fragmentation model [18]). 

According to the consideration mentioned in sec. 2, if the spin of the charmed hadron 

is same as the spin of the charm quark produced in the subprocess, the model (A) 

might be a reasonable scenario " 5 . Concerning the model (B), we follow the analysis 

for A production [18] and apply it to A+ / D* + productions. 

Since the charm quark is heavy, the Bjorken-x of the parton taking part in 

charm quark pair production in the subprocess is not very large for pp collisions with 

an appropriate energy. In this charm quark is expected to be dominantly 

produced via gluon-gluon fusion because of the rather large gluon distribution. 

However, since valence quark densities become larger in large x region, we could 

not neglect a contribution from a quark-anti-quark annihilation qq — > cc in high 

" 5 From HERMES experiment, even the case of A baryon which include lighter s quark, a frag- 
mentation function based on the naive-quark-parton model seems to be reasonable scenario [19]. 
Therefore, the model (A) must not be also an unreasonable scenario for A+/_D*+ productions. 




(11) 



p T region. Therefore, in order to get a good information on the gluon polarization, 
we must find the kinematical region where gg — > cc dominates over qq — > cc. To do 
so, we calculated the p T distribution of the cross sections dAa/dp T , da/dp T and 
spin correlation asymmetries D LL for D* + and A+ productions at v /s=200GeV 
( v /s=500GeV) as shown in Fig. 2 (Fig. 3) and Fig. 4 (Fig. 5), respectively. In the 
left panel of these figures, we present the results only for the model (A), since the 
purpose of this analysis is to compare a contribution of gg — > cc with the one of 
qq — > cc. Since in the p T region of 3 < p T < 5 GeV, the cross section, dAa/dp T 
and da/dp T , of gg — > cc is much larger than the one of qq — > cc, the _D LL for the 
sum of gg — > cc and gg — ■> cc do not change much from the one for gg — > cc alone. 
Therefore, we can say that a contribution from the subprocess gg — > cc can be 
safely neglected in the region of 3 < p T < 5 GeV. 

We also calculated D LL as a function of rj for 3 kinematical regions of rj, i.e. (i) 
—0.5 < rj < 0, (ii) —0.5 < rj < 0.5, and (iii) < rj < 0.5. However, since the result 
for the region (i) distinguishes most clearly a contribution of gg — > cc from the one 
of gg — > cc, we show only the result for the region (i) in Fig. 2 ~ Fig. 5. For this 
limited p T region (3 < p T < 5 GeV), D LL for D* + and A+ production are shown 
in Fig. 6 ~ Fig. 7, respectively, at y/s = 200GeV and v /i=500GeV, where we see 
a rather big model-dependence of polarized gluons and also the spin-dependent 
fragmentation functions. 

Furthermore, to know how Dll is sensitive to the behavior of polarized gluons 
in the proton, we calculated the statistical sensitivities of D LL , i.e. SD LL , which is 
defined by the following formula [20]; 



To numerically estimate the value of SD LL , we used the following parameters: the 

beam polarization; P =70%, a integrated luminosity; L=320pb~ 1 (800pb~ 1 ) for 

y/s = 200 (500) GeV [3], the trigger efficiency; e = 10% for detecting produced 

charmed hadron events, decay asymmetry parameter; a = 1.0 (—0.98) for D* + 

(A+) decay and a branching ratio; b D *+ = Br(D* + -> D°7i + )Br(D° -> R-ti + ) ~ 

2.5 x 10~ 2 , b A + = Br(A+ -> A7r+)Br(A -> pir') ~ 5.8 x 10~ 3 [10] » 6 . a denotes 

" 6 In order to determine the polarization of D* + and A+, observation of the angular distribution 
for decay channels of these particles is necessary. Practically this could be done by detecting 
charged particles in final the state. 




(12) 



the unpolarized cross section integrated over a corresponding r\ region. The factor 
3 in Eq. (12) is an acceptance factor for our processes [20]. In Figs. 6 and 7, 
statistical sensitivities, 5D LL , are attached only to the dashed line of D LL which 
were calculated using the GRSV01 parametrization model of polarized partons and 
the non-relativistic fragmentation model (model (A)) * 7 . As shown here, statistical 
sensitivities 5Dll are so small that the processes must be feasible for measuring 
the D LL . 

From these results, we see that the 77 distributions of D LL are effective for test- 
ing the model of not only polarized gluon distributions but also spin-dependent 
fragmentation functions for all cases (D* + production and A+ production) calcu- 
lated at v/i =200 GeV and 500 GeV. Especially, rj distributions of D LL for D* + 
productions are quite promising, though the magnitude of D LL becomes smaller 
with increasing center of mass energy from y/s =200 GeV to 500 GeV. 
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Figure 2: The cross sections (left panel) and Dll (right panel) for D* + productions 
as a function of p T at i/i=200 GeV. Here, gg, qq and gg + qq shown in the paren- 
thesis at explanatory notes correspond to the subprocess taken into consideration. 



5. Conclusion 

In order to get information about polarized gluons in the proton, we proposed 
two charmed hadron production processes, i.e. polarized D* + meson productions 
and polarized A+ baryon productions, which will be observed in the forthcom- 
ing RHIC experiments. As described in Introduction, the processes contain two 

" 7 Note that as shown in Eq.(12), SDll does not depend on both the model of polarized gluons 
and the model of fragmentation functions. 




Figure 3: The same as in Fig. 2, but for v /i=500GeV. 




Figure 4: The same as in Fig. 2, but for A+ production. 



production mechanisms, gluon-gluon fusion (gg — > cc) and quark-anti-quark an- 
nihilation (qq — > cc). Thus, to study the gluon polarization in the proton, it is 
necessary to find the kinematical region where gg — > cc dominates over qq — > cc. 

From the numerical calculation at the lowest order of QCD, we found that the 
i] distribution of Dll in the limited p T region (3 < p T < 5 GeV) is quite promising 
for testing not only the model of polarized gluons in the proton but also the model 
of spin-dependent fragmentation functions. However, it should be noted that the 
assumption on the spin-dependent fragmentation function might be somewhat too 
simple in the present analysis. In order to obtain more reliable prediction, further 
investigation on spin-dependent fragmentation functions is necessary. For the A 
baryon production, similar analysis was performed by D. de Florian et. al [21]. Our 
analysis might be complementary to their analysis but more effective for extracting 
information of the polarized gluon in the proton, since the non-relativistic quark 
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Figure 5: The same as in Fig. 3, but for A+ production. 
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Figure 6: The Dll for D* + productions as a function of i] at ^=200 GeV(left 
panel) and ^=500 GeV(right panel) 



model works better for A+ (D* + ) than A and furthermore the separation of gluon 
fusion and qq annihilation is easier for our processes. 

Though the present calculation is confined to the leading order, the results are 
interesting and we hope that our analysis could be tested in the forthcoming RHIC 
experiments. 
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